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ABSTRACT 


The  purpose  of  this  Investigation  was  to  determine  analytically  the 
molecular  incidence  rates  on  the  various  surfaces  of  the  space  simulator. 

The  mathematical  model  of  the  space  simulator  consists  of  two  concentrically 
located  spherical  surfaces  corresponding  to  the  chamber  wall  outside  and  to 
the  vehicle  inside.  The  analysis  is  divided  into  two  phases. 

Phase  1  deals  with  conditions  of  uniform  pumping  and  outgassing  on 
each  surface.  A  square  pulse  of  outgassing  is  assumed  to  occur  and  the  total 
number  of  molecular  hits  is  determined  for  each  surface.  These  results  are 
extended  to  a  steady  state  uniform  flow  and  the  molecular  incidence  rates  are 
obtained. 

Phase  11  is  an  analysis  of  the  nonuniform  case  introduced  by  surface 
discontinuities  of  axiosymmetric  type.  The  procedure  followed  is  to  deter- 
mine  the  various  probabilities  that  a  molecule  leaving  a  surface  will  hit 
spherical  zones  on  either  surface.  Using  these  probabilities  the  molecular 
incidence  rates  can  be  determined.  While  the  probabilities  have  been  found, 
analytic  expressions  for  molecular  incidence  determination  are  very  difficult 
to  obtain. 

The  problem,  therefore,  was  approached  from  a  computer  point  of  view 
enabling  us  to  obtain  computer  solutions  with  a  high  degree  of  accuracy. 
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SYMBOLS 


Subscripts  v  and  c  refer  to  vehicle  and  chamber  wall,  respectively.  In  case 
we  specify  a  portion  on  either  of  these  two  surfaces,  v  and  c  will  be  employed 
accordingly. 

a^  -  pumping  coefficient  on  surface  x  (where  x  «  v  or  c) 

P  »  probability  that  a  molecule  leaving  surface  x  will  hit  surface  y 

(«cor  v)  without  an  intermediate  collision.  Clearly  P  *  1,  P  *0 

H  ■total  number  of  molecular  hits  against  the  non- pumping  surface  y  due 

to  a  flow  originated  at  surface  x. 

H  ■  total  number  of  molecular  hits  against  the  non-pumping  portion  of 

surface  x  (x  ■  v  or  c) 

M  (n)  ■  number  of  molecules  pumped  out  on  surface  x  during  the  interval 

*  nT<t<(tH-l)T,  n«l,2,.  .  • 

■  total  number  of  molecules  pumped  out  on  surface  x 

f  (n)  ■  number  of  molecules  leaving  surface  x  during  the  interval 

*  nT<t<(n4-l)T ,  n“l ,  2 ,  •  •  • 

F  ■  steady  state  luiiform  flow  per  unit  area  per  unit  time  originating  at 

surface  x 

■  rate  of  molecular  incidence  on  surface  x  per  unit  area 

r  ■  radius  of  surface  x 

X 

y  ■  angle  between  the  tangent  vector  to  the  vehicle  and  the  radial  vector 

both  vectors  drawn  from  a  point  on  the  chamber  wall. 

In  considering  the  probability  P  of  a  point  p  (see  Figures  2, 3, 4, 5, 6, 7) 

xy 

6^  ■  latitude  of  the  point  p 

■  boundary  of  the  surface  region  on  y  illxaminated  from  p 
P 

6  ■  latitude  on  surface  y 

o 

■  subregion  of  the  illtuninated  region  on  y  bounded  by  and/or  6^ 

o  for  which  the  latitude  Q  of  any  point  in  that  subregion  satisfies 

the  relation  &>Q 
—  o 


lx 


AEDC<TDR-63-88 


xy'  i’ 


p(xPn 


Sphtrical  Zoii«  Probobi 

«  probability  that  a  molecule  departing  from  a  point  on  surface  x 
with  latitude  6.  will  hit  the  surface  region  on  surface  y  bounded 
below  by  the  laHtude  0..  This  region  may  or  may  not  be  illuminated 
completely  from  that  point  (see  Figures  8  -  15). 

-  probability  that  a  molecule  departing  from  a  spherical  zone 
will  hit  the  spherical  zone  Yj. 

-  number  of  molecules  leaving  the  spherical  zone  x.  during 

nT<t<(rH"l)T 

=  i*th  spherical  zone  on  surface  X 

=  spherical  zone  on  the  chamber  wall  whose  upper  boundary  is  the 
discontinuous  line. 


Single  Surface  Discontinuity 


=  a  for  all  i 

V 


=  a  l<i<N- 


=  a  n,+l<i<2N 
C2  (^1  -- 

=  flow  rate  per  unit  area  introduced  on  the  part  of  the  chamber 
wall  for  which  a  »  a 


X,  Y,  Z 
r,  e,  $ 


=  flow  rate  per  unit  area  introduced  on  the  part  of  the  chamber 
wall  for  which  a  =  a 

Ci  c^ 

flow  per  unit  area  introduced  on  the  vehicle 
=  usual  cartesian  coordinate  system 

-  usual  spherical  coordinate  system 
=  complement  of  6 

=  same  as  $ 

-  spherical  zone  on  either  surface  bounded  by  the  latitudes  6. 
and  ©j 

=  latitude  of  discontinuous  line  on  chamber  wall 


N  »  number  of  spherical  zones  of  a  hemisphere 

N2  =  7/^ 

$1  *  lower  latitude  of  the  I*th  spherical  zone 

Liin  ■  maximum  latitude  of 

Lim  »  lower  latitude  of 


X 
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1.0  INTRODUCTION 


As  a  necessary  consequence  of  the  space  age.  It  has  become  mandatory 
to  plan  for  ground  testing  of  space  vehicles,  components  and  materials.  The 
duplication  of  all  the  spatial  environments  is  obviously  impossible.  Hany 
papers  Jiave  been  written  defining  those  environments  which  could  be  simulated 
and  the  degree  to  which  simulation  is  required.  Among  those  listed  as  neces¬ 
sary  and  feasible  is  the  low  pressure  environment  of  outer  space.  It  is 
frequently  stated  that  a  pressure  of  10*^  Torr  is  sufficient  to  simulate  the 
thermal  radiation  environment  and  the  high  electrical  impedance  of  space. 

This  statement  alone  does  not  define  where  the  pressure  is  measured,  nor  in 
fact  what  is  meant  by  pressure  in  the  highly  directional  conditions  of  both 
true  space  and  a  space  environmental  simulator.  A  serious  evaluation  of  the 
nature  of  the  molecular  flux  in  both  cases  is  necessary  in  order  to  evaluate 
what  is  meant  by  simulation. 


1.1  THE  VEHICLE  IN  SPACE 


Two  altogether  different  gaseous  sources  surround  the  vehicle  in  space. 
First,  there  is  the  natural  environment  which  was  there  prior  to  the  entry 
of  the  space  vehicle,  and  secondly,  the  gas  sources  which  have  been  contribu¬ 
ted  by  the  vehicle. 

The  atmosphere  of  the  planet  earth  at  altitudes  above  1500  km  consists 
primarily  of  hydrogen  ions,  decreasing  in  density  with  Increase  in  altitude. 
Below  this  level,  the  atmosphere  consists  principally  of  oxygen  and  nitrogen 
in  atomic,  ionized,  or  molecular  states,  depending  on  the  altitude.  Between 
200  and  1500  km  the  density  decreases  with  increasing  altitude  from  about 
10^0  to  10^  particles  per  cubic  centimeter.  The  corresponding  mean  free  paths 
range  from  a  few  miles  to  several  hundred  thousand  miles,  distances  which  are 
huge  relative  to  the  dimensions  of  any  space  craft.  The  velocity  of  these 
molecules  ranges  from  1  to  2  km/ sec.  Yet,  the  velocity  of  an  earth  satellite 
ranges  from  7.5  to  11  km/ sec.  The  vehicle  then  appears  to  be  rushing  through 
a  relatively  stationary  low  density  gas  medium,  most  of  the  collisions  being 
on  the  leading  edge  of  the  space  craft.  We  note  then  that  the  molecular  in¬ 
cidence  rate,  and  the  pressure,  depends  on  the  velocity  of  the  vehicle  and 
varies  over  the  surface  of  the  vehicle. 


In  addition  to  the  planet  atmosphere,  the  sun  is  thought  to  emit  a 
variable  stream  of  protons  and  electrons  (solar  wind)  having  a  density  of 
600  to  10,000  protons  per  cubic  centimeter  near  the  earth.  The  velocity  of 
these  protons  varies  from  500  to  2Q00  km/ sec.  depending  on  the  solar  activity. 
The  incidence  rate  is  from  5  x  10^^  to  2  x  10^^  particles  per  second  per 
square  centimeter,  which  is  comparable  to  the  Incidence  rate  of  molecules 
from  the  planet  atmosphere  at  high  altitudes.  Since  the  velocity  of  the 
solar  wind  is  so  much  higher  than  that  of  the  vehicle,  the  collisions  are  all 
on  the  side  facing  the  sun.  The  extremely  high  energy  is  sufficient  to  cause 
sputtering  or  other  surface  phenomena.  Again  we  note  that  there  is  no  uni¬ 
formity  over  the  surface  of  the  vehicle,  nor  can  we  simply  define  a  pressure 
equivalent  of  altitude. 
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In  addition  to  the  environmental  gases,  the  vehicle  itself  is  a  major 
source  of  gas,  in  the  form  of  outgassing  and  leakage.  This  gas  is  given 
off  by  the  vehicle  with  thermal  energies  relative  to  the  vehicle.  The 
outgassing  occurs  in  all  directions,  but  may  well  vary  over  the  surface  of 
the  vehicle,  depending  on  the  variations  in  the  vehicle  structure  and  mater¬ 
ial  composition.  The  outgassing  products  will  travel  many  miles  in  a 
straight  line  before  making  collisions  with  an  environmental  molecule,  hence 
for  all  practical  purposes,  the  molecules  which  leave  the  vehicle  are  per¬ 
manently  gone.  Space  may  therefore  be  considered  as  pumping  molecules  as 
fast  as  they  are  produced.  This  pumping  action  is  not  related  to  the  con¬ 
ventional  concept  of  effusion  velocity  of  an  orifice  of  a  fixed  speed  per 
unit  area,  as  no  fixed  area  of  space  can  be  defined  as  the  pumping  surface. 

In  almost  all  space  vehicle  situations,  the  outgassing  of  the  vehicle 
represents  a  far  greater  density  than  that  contributed  by  either  the  planet 
environment  or  the  solar  wind.  Hence,  the  gas  density  surrounding  the 
vehicle  is  almost  entirely  made  up  of  the  composition  of  the  outgassing  pro¬ 
ducts.  Yet  this  gas  is  all  moving  away  from  the  vehicle  and  does  not  con¬ 
tribute  to  the  gas  returning  to  the  vehicle.  The  returning  gas  is  still 
defined  by  the  velocity  of  the  vehicle  and  the  normal  environment  through 
which  it  is  moving. 

1.2  the  vehicle  in  a  chamber 

When  a  vehicle  is  placed  in  a  space  simulation  chamber,  it  is  again 
exposed  to  two  different  sources  of  gas:  1)  The  natural  environment  or 
background  of  the  chamber,  and  2)  The  gas  contributions  due  to  the  out¬ 
gassing  and  leakage  of  the  vehicle. 

For  true  duplication  of  the  atmosphere  of  space,  the  chamber  back¬ 
ground  should  match  that  of  the  upper  atmosphere  in  gas  type,  density  and 
energy.  In  addition,  it  should  be  properly  oriented  to  account  for  the 
relative  velocity  of  the  vehicle  through  the  satationary  gas,  and  the  solar 
wind.  While  such  conditions  are  not  completely  impossible,  they  are  suffi¬ 
ciently  difficult  and  expensive,  that  realistic  evaluations  are  required  to 
determine  the  necessity  of  such  duplication.  In  the  normal  space  chamber, 
the  residual  gas  is  made  up  of  the  outgassing  products  of  the  walls  of  the 
chamber  and  leaks  which  may  exist  to  the  outside  atmosphere  or  to  various 
internal  fluid  sections  such  as  cryogenic  lines.  The  composition  of  this 
gas  will  vary  in  different  chambers,  and  may  include  such  gas  types  as: 
water  vapor,  hydrogen,  carbon  dioxide,  carbon  monoxide,  nitrogen  and  helium. 
The  gas  is  not  ionized,  and  in  general  is  at  a  temperature  of  about  lOO^K. 

In  addition,  there  may  be  directional  effects  in  the  event  that  discontinu¬ 
ities  exist  on  the  chamber  wall.  For  example,  the  presence  of  a  large  area 
solar  simulator  may  contribute  a  high  outgassing  rate  from  a  particular 
section  of  the  chamber  wall. 

When  the  vehicle  is  in  a  chamber,  just  as  when  it  is  in  space,  the 
outgassing  of  the  vehicle  is  the  predominate  gas  present.  However,  in  the 
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case  of  the  chamber  an  altogether  different  situation  exists.  In  order  to 
duplicate  the  effectively  infinite  pumping  speed  of  outer  space  the  chamber 
must  have  walls  which  are  perfect  condensers  for  all  incident  molecules. 

Since  such  perfect  condensers  are  not  presently  feasible »  a  great  majority 
of  the  outgassing  from  the  vehicle  is  reflected  from  the  chamber  wall  and 
returns  to  the  vehicle.  In  almost  all  Instances,  this  returning  gas  flux 
is  far  larger  than  the  original  chamber  background.  Again  we  note  that 
this  gas  flux  does  not  have  the  characteristic  properties  of  the  permanent 
space  background.  The  gas  is  predominately  water  and  nitrogen.  The  tempera¬ 
ture  is  between  100  and  300^K.  Anomalous  directional  effects  may  exist  if 
large  areas  of  non  pumping  walls  are  present. 

The  foregoing  differences  between  a  chamber  and  space  may  not  be  of 
importance  in  the  superficial  condition  of  simulating  for  heat  transfer  and 
electrical  discharge  effects.  However,  if  any  surface  phenomena  is  involved, 
such  as  friction  or  sputtering,  entirely  different  results  may  be  obtained 
in  the  two  cases. 

1.3  MOLECULAR  KINETICS 

Considering  the  foregoing  differences  between  true  space  and  so-called 
space  simulators,  it  is  apparent  the  simple  question  "What  altitude  has  been 
simulated?"  is  almost  impossible  to  define.  In  the  simplest  case  it  creates 
a  serious  controversy  of  philosophy,  measuring  equipment,  and  people.  The 
Molecular  Kinetics  program  was  started  to  study  the  various  problems  of 
gas- vehicle  collisions  in  space  and  in  a  chamber,  to  produce  a  yardstick 
for  measuring  simulation,  and  to  define  the  effects  or  limitations  of  vari^ 
ous  chamber  parameters  on  the  degree  of  simulation  achieved.  This  program 
is  divided  into  a  number  of  different  phases,  each  covered  in  a  different 
report.  The  common  denominator  of  all  phases  of  the  Molecular  Kinetics 
program  is  the  dealing  with  highly  directional  gas  flows,  in  terms  of  measur¬ 
ing,  controlling,  or  defining  simulation. 

In  order  to  derive  a  language  to  be  used  in  the  interpretation  of 
space  simulation,  and  which  is  compatible  with  the  parameters  present  in 
space  simulator,  it  is  first  necessary  to  define  the  ground  rules  to  be 
used.  We  first  note  that  the  vehicle  is  the  point  of  interest  for  dis¬ 
cussing  simulation,  and  further  it  is  the  molecules  arriving  at  the  vehicle 
which  are  of  importance.  Any  measurement  or  consideration  of  molecular 
fluxes  or  pressures  on  the  chamber  wall  are  of  importance  only  in  so  far  as 
they  can  be  interpreted  at  the  vehicle  proper.  The  second  important  factor 
to  be  determined  is  the  term  to  be  used  for  discussing  gas  particles  in  the 
simulator.  Three  different  terms  are  immediately  apparent:  pressure, 
density,  and  particle  flux. 

While  density  is  very  meaningful  for  the  random  gas  in  space  in  the 
absence  of  a  vehicle,  it  has  little  meaning  for  the  vehicle  when  a  large 
amount  of  outgassing  products  from  the  vehicle  are  included  in  the  measure¬ 
ment.  The  density  concept  does  not  define  the  direction  of  the  gas  and 
hence  cannot  be  related  to  the  molecules  striking  the  vehicle,  either  in 
space  or  in  a  chamber. 
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A  consideration  of  Che  earlier  discussions  on  the  extreme  directional 
effects  both  in  space  and  in  a  chamber  indicates  that  molecular  incidence 
and  energy  level  are  the  important  factors  in  assessing  simulation.  It  is 
apparent  that  there  is  a  considerable  difference  in  the  temperatures  in 
space  and  in  a  simulator.  A  theoretical  assessment  of  the  temperature  of 
the  molecules  received  by  the  vehicle  in  a  space  chamber  is  almost  impossi* 
ble  due  to  the  variations  in  the  temperature  within  the  chamber  and  the 
differences  in  the  accommodation  coefficient  for  different  gases.  While 
pressure  is  related  to  both  energy  and  incidence  rate,  the  true  measurement 
of  pressure  is  likewise  a  very  difficult  problem.  From  both  the  analysis 
and  the  measurements  point  of  view  it  is  far  easier  to  consider  the  rate  at 
which  molecules  are  striking  a  surface.  Except  for  cases  of  extreme  energy 
(such  as  the  solar  wind)  the  surface  reactions  are  more  dependent  on  the 
gas  type  than  on  the  variations  of  energy  which  exist.  For  these,  and 
other  reasons  which  will  be  more  apparent  in  the  different  reports,  the 
concept  of  molecular  flux  has  been  chosen  for  all  analysis.  Molecular  flux 
is  defined  as  the  number  of  molecules  per  square  centimeter  per  second 
arriving  or  leaving  a  given  surface. 

1.4  MOLECULAR  FLUX  DISTRIBUTION 

This  part  of  the  Molecular  Kinetics  program  has  to  do  with  the  mathe¬ 
matical  analysis  of  computing  the  molecular  flux  on  every  part  of  the  space 
vehicle  and  the  chamber  walls  as  a  result  of  outgassing  occurring  from 
either  the  vehicle  or  the  chamber  wall.  The  program  further  considers  the 
effects  of  a  surface  discontinuity,  such  as  a  solar  region,  on  the  flux 
distribution.  In  solving  this  problem  probability  equations  were  first  der¬ 
ived  for  an  axial  symmetry,  and  equations  were  generated  for  the  limiting 
conditions  of  every  possible  configuration.  In  all  cases  a  cosine  distribu¬ 
tion  was  assumed  for  the  leaving  molecules.  The  general  problem  of  a  single 
surface  discontinuity  was  then  handled  by  means  of  a  computer  program.  The 
program  considers  the  case  of  a  spherical  vehicle  in  a  spherical  chamber, 
one  of  which  contains  a  spherical  cap  of  different  properties.  The  surfaces 
are  then  divided  into  circular  sections  about  a  ccxnmon  axis  congruent  with 
the  axis  of  the  single  surface  discontinuity.  A  value  of  sticking  coeffi¬ 
cient  and  outgassing  rate  is  assigned  to  each  surface.  The  computer  program 
counts  the  molecules  arriving  and  leaving  each  and  every  surface  until  a 
steady  state  flow  condition  is  obtained*  This  final  value  then  defines  the 
molecular  flux  to  the  different  parts  of  the  system  as  a  result  of  the 
assumed  outgassing  rates. 
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2.0  COSINE  DISTRIBUTION 


In  the  molecular  migration  from  one  surface  to  another,  it  will  be 
assumed  that  all  molecules  follow  a  "Cosine  Distribution";  that  is,  the 
probability  that  a  molecule  will  travel  in  a  certain  direction  when  it  leaves 
from  a  point  on  a  surface  is  proportional  to  the  cosine  of  the  angle  between 
this  direction  and  the  normal  to  the  surface  at  the  point  in  question.  Thus, 
molecules  have  the  greatest  tendevtcy  to  leave  a  surface  in  a  direction  normal 
to  the  surface.  The  statistical  concept  of  this  law  is  presented  in  the 
following  discussion. 

Let  dA.  be  a  surface  element  on  a  surface  and  let  us  assume  that  it 
is  required  to  find  the  probability  that  the  molecular  flow  emitted  from  that 
surface  element  will  hit  a  specified  area  A2  on  a  surface  S2.  If  the  element 
area  dA^  is  taken  sufficiently  small  we  may  speak,  in  a  sense,  of  a  point  p 
instead  of  the  element  dA.  as  if  the  total  flow  from  dA,  was  emitted  from  the 
point  p.  If  a  sphere  S  of  arbitrary  radius  r  (unit  radius  will  be  convenient) 
is  constructed  with  center  the  point  p,  the  directional  vector  from  the  point 
p  to  the  boundary  of  area  A*  will  describe  a  closed  curve  s  on  the  surface  S. 
Let  A  be  the  area  enclosed  by  s.  If  dA  is  an  area  element  of  A  we  have 

2 

dA  =  r  cosed$d0 


Let  6  =  f(<I>)  be  the  equation  of  the  curve  s.  The  angle  between  the 
normal  to  at  p  and  the  vector  to  any  point  of  s  from  p  is  Since  we 

assume  a  cosine  distribution  it  follows  that  the  molecular  flow  through  dA 
will  be 


dA  Kcos(^-6)  =  Kr^cos0sin0d<l>d0 


-  Kr  sinoecosoed^da 


where 

Ot  -  K  is  a  proportionally  constant 


Thus 

(Total  Flow  thru  A) 


f-f(P) 

2 

/d^/  Kr  cosOsinUda 


On  the  other  hand 

(Total  Flow  thru  hemisphere) 


=  /d3/ 


It 

*2 

2  2 
Kr  cosOtsinCifda  »  Kitr 


o 


5 


4 


AEDC-TDR.63*88 


Then  the  probability  P  of  hitting  the  area  A  from  the  point  p  is 


p  ^  (Total  flow  thru  A) _ 

(Total  flow  thru  hemisphere 


f-fO) 

1  ^ 

—  /d3/  cosOCsivCXdcx 

^  n 


Thus,  the  probability  that  a  molecule  emitted  from  a  point  on  a  sphere  of 

radius  r  will  hit  an  inside  sphere  of  radius  r  is 
c  V 


1  ^  2 

p  a  “  J  dp/  cosOesinOda  «  sin  y 
R  0  o 
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3.0  UNIFORM  PUMPING 


The  objective  of  this  section  is  to  determine  the  rate  of  molecular 
incidence  on  either  surface  under  the  most  simplified  conditions.  In  this 
analysis,  the  most  simplified  model  of  space  simulator  is  considered;  i.e«, 
a  spherical  test  vehicle  located  concentrically  in  a  spherical  test  chamber. 
The  walls  of  both  the  vehicle  and  the  chamber  are  assumed  to  have  uniform 
properties  and  the  pumping  coefficients  on  each  of  these  two  surfaces  are 
assumed  uniformly  distributed.  To  be  more  specific,  if  we  define  the 
pimping  coefficient  a^  by: 

a  =  A,/A 

X  1 

where  A.  is  the  pumping  portion  of  area  A,  a  remains  constant  throughout 
the  surface  x.  Furthermore,  in  order  to  covir  every  special  case  we  consider 
the  general  case  where  ^  ^  0. 

The  analysis  is  made  on  the  basis  of  a  pulse  of  molecules  being  admitted 
uniformly  from  either  the  vehicle  or  the  chamber  wall  or  both  surfaces.  The 
results,  however,  are  equally  applicable  to  a  steady  state  flow  as  it  will  be 
shown  in  a  derivation  below. 

With  this  mathematical  model  we  have 
P  -  sin^y  -  P,  P  -  1 

3.1  DIRECTIONAL  FLOW  (PULSE) 

3.1.1  Differance  Equations 

The  following  difference  equations  are  easily  derived 


-  P(l-a^)f^(n)  (3.1) 

f  (iH-2)  -  (1-a  )(1-P)f  (nfl)+(l-a  )f  (nfl)  (3.2) 

C  C  C  C  V 

Substituting  3.1  into  3.2  we  obtain 

f^(iri-2)-(l-a^)(l-P)f^(iri-l)-(l-a^)(l-a^)Pf^(n)  -  0  (3.3) 

The  general  solution  of  3.3  Is 
f^(n)  -  Cj^r^"+C2’^2" 


where  C.  and  C2  are  arbitrary  constants  which  have  to  be  defined  from  the 
boundary  conditions,  and  r^  and  r^  are  the  roots  of  the  characteristic  equation 
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t(r)  -  r^-(l-a^)(l-P)r-(l-a^)(l-a^)P  -  0 


(3.5) 


Boundory  Conditions 


If  the  flow  is  Introduced  on  both  surfaces,  f  (I)  will  not  depend  on 
f  (0)  exclusively  but  on  both  f  (0)  and  f  (0).  If  ^he  flow  is  introduced 
on  surface  x  only  then  f  (0)»0  for  y  9^  x.  In  any  case  we  assume  f  (0)  and 
f^(l)  given  covering  thu^  all  special  cases.  ^ 

Then  from  3.4  by  setting  n  »  0,1  we  obtain  the  following  system  of 
equations 


f^(0)  » 

f  (U  -  C^r  +  Cjrj 


(3.6) 


Determining  the  constants  and  from  3.6  and  putting  them  in  4) 
we  obtain  finally 


fc(n) 


£c(0)r^r2(ri"-^-r2---S 

^2~^L 


fpCDCri^-ra") 

r2-ri 


(3.7) 


However  the  number  of  molecules  leaving  the  chamber  wall  during 
nT<t<(n+l)T  Is  the  same  as  the  number  of  molecular  hits  on  the  nonpumping 
part  of  the  chamber  wall  during  nT<t<(nfl)T.  Hence,  the  total  number  of 
molecular  hits  on  the  chamber  wall  Is 


H  =  E  f-(n) 
c  n=l 


(1-rj^)  (l-r2) 


(3.8) 


Similarly,  from  equation  3.1  and  3.8  we  obtain 


»v  ■  fc<”) 


.  P(l-a  ) 

V 

-  P(l-a  ) 

V 


fc(l)-rir2fc(0) 

(l-rj^)(l-r2) 


+  f^^o) 


fc(0)[l-(ri+r2)]  +f^(l) 

(l-r^)(l-r2) 


-  P(l-a  ) 

V 


(3.9) 
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Another  form  of  equations  3.8  and  3.9  in  terms  of  f  (0)  and  f  (0)  can 
be  obtained  by  making  the  following  substitution  ^  ^ 

f^(l)  «  (l-a^)(l-P)f^(0)+(l-a^)yQ^ 

-  (ri+r2)f^(0)+(l-a^)f^(0) 


Then  3.8  and  3.9  become 

(r^+r2-r^r2)f^(0)+(l-a^)f^(0) 
“  (1-r^)  (l-r2) 

f  (0)+(l-a  )f  (0)  . 

V  (l-rj^)(l-r2) 


(3.10) 

(3.11) 


3.1.2  Special  Casts 


Two  special  cases  are  of  Interest,  namely, 

f^(l)  -  (l-P)(l-a^)f^(0)  =  (r^+r2)f^(0),  f^(0)  -  0 

f^(l)  -  (l-a^)f^(O)  ,  f^(0)  -  0 


Substituting  these  two  equations  in  equations  3.8  and  3.9  or  Just  simply 
setting  f  (0)  or  f  (0)  equal  to  zero  In  equations  3. 10  and  3.1]l  we  obtain  the 
following  special  cases. 


H 


ri+rz-rirz 


P(l-a  ) 

-I  V 

cv  ■“  (l-rj^)(l-r,)  '■c 

(l-Hc) 


~^1^2 


£e«» 

(3.12) 

y  £,(o)  -  0 

l-h 

o 

o 

(3.13) 

< 

o 

(3.14) 

y  £,.(0) .  0 

£vC0)  ) 

(3.15) 
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3.1.3  Moltcular  Pumping 

The  number  of  molecules  pumped  on  either  surface  during  nT<t<(n4-l)T  is 
given  by: 


a 


M  (n)  =  7-^ 
c'  1-a 

c 

(3.16) 

a 

f  (n) 
v' 

(3.17) 

V 


The  total  number  of  molecules  pumped  is 


M  =  I  M  (n) 
n=l  c 


H 


c 


uu 

M  =2:  M  (n) 
V  n=l  V 


H 


V 


(3.18) 


(3.19) 


Note  that  we  must  have 

M  +  M  =  f  (0)  +  f  (0) 

C  V  c  V 

since  eventually  all  the  introduced  molecules  must  be  pumped  out.  Tliis  can 
be  verified  to  be  the  case, 

3.2  STEADY  STATE  UNIFORM  FLOW 

The  results  obtained  in  3.1  can  be  extended  to  the  case  where  a  steady 
state  molecular  flow  is  established  instead  of  a  pulse  under  the  following 
assumptions. 

a)  The  rate  of  molecular  incidence  is  determined  after  a  long 
period  of  time  since  the  time  when  the  uniform  steady  state  flow  was  initiated. 

b)  Time  duration  between  two  successive  hits  of  each  molecule  is 
taken  to  be  a  statistical  average  time  interval  t. 

Consider  now  a  square  pulse  of  flow  into  the  system  of  time  duration 
starting  at  t  «  0  and  suppose  that  we  want  to  determine  the  number  of  molecular 
hits  during  the  time  interval  nT<t<(nfl)T. 

According  to  assumption  b) ,  the  number  of  molecule  hits  during  the  time 
interval  nT<t<(rtfl)T  will  depend  on  how  many  molecules  are  present  in  the 
space  simulator  during  the  time  interval  (n-l)T<t<nT.  Thus,  our  problem  is 
to  determine  the  number  of  molecular  hits  f(n)  during  nT<t<(n+l)T  due  to  a 
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square  pulse  established  during  the  time  interval  0<t<T|  or  saying  the  same 
thing  in  a  different  way,  to  find  the  niunber  of  molecular  hits  during  0<t<T 
due  to  aquare  pulse  established  -nr  ago.  If  we  had  two  identical  pulses 
established  at  t  ■  nx  and  t  ■  -(nfm)T,  each  pulse  of  duration  t  then  by  super¬ 
position,  the  total  number  of  molecular  hits  during  ()<t<T  due  to  these  two 
pulses  will  be  f(n)+f(nfm).  Now  we  are  ready  to  pass  to  a  steady  uniform 
flow.  If  the  flow  is  considered  to  consist  of  square  pulses,  the  total  number 
of  molecules  present  during  the  time  interval  (KtKj  is 
00 

Z  f(n)  where  N  is  very  large 
n=0 


If,  however,  f(n)  converges  very  rapidly  (as  indeed  is  the  case)  we  may 
replace  N  by  infinity.  Then 

N  00 

£  f(n)^£  f(n) 
n=0  n=0 

Incidentally,  it  is  obvious  that  any  transient  conditions  of  the  flow  at  the 
beginning  will  have  negligible  effects. 

00 

Now  let  H  be  the  number  of  molecular  hits  during  0<t<T  due  to 
molecules  present  during  -T<t<0.  Then  the  rate  of  molecular 
incidence  during  0<t<T  is  H. 


3.2,1  Rate  Determination  of  Molecular  Incidence 
Clearly 

f  (0)  =  (l-a^)F  T4iir  ^ 
c  c  c  c 

f  (0)  =  (l-a  )F  T4jtr  ^ 

V  V  V  V 

Substituting  these  in  equations  3.10  and  3.11  we  obtain 

_  (ri+r2-r^r2)F^+(l-a^)F^P 
\  “  (l-rj^)(l-r2) 


(3.2.1) 


(3.2.2) 
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3.2.2  Sptciol  Casts  of  Rato  Dotormi nation 


a) 


»  0  (flow  originated  at  chamber  wall) 


(r,+r«-r.r^)F^ 
R  ^  1  ^  ^2 _ c 

^  (l-rpCl-r^) 


R  == 


(1-a  )F 
'  c  c 


V  (l-r^d-r^) 


If  a^  =  0,  equations  3.2.3  and  3.2.4  reduce  to 


R  =  R  = 
c  V 


(1-a  )F 
^  c  c 

a 


b) 


=  0  (flow  originated  at  vehicle) 


(1-a  )PF 

^  V  V 

C  “  (l-rj^)(l-r2) 


R  = 


R  = 


ri2.F 

1  2  V 


V  “  (l-r-)(l-r2 


(3.2.3) 

(3.2.4) 


(3.2.5) 


(3.2.6) 

(3.2.7) 


If  a^  =  0,  equations  3.2,6  and  3.2.7  reduce  to 
PF 

R  =  — ~ 
c  a 

c 

(1-a  )PF 

«  _  _ _ y 


(3.2.8) 


(3.2.9) 


3.3  RESUME 

The  results  obtained  by  admitting  a  square  pulse  of  molecular  flow  on 
either  surface  are  represented  by  equations  3.8  and  3.9  or  3.10  and  3.11  and 
special  cases  are  represented  by  equations  3.12  -  3.15.  The  total  number  of 
molecules  pumped  followed  Immediately  from  equations  3,18  and  3.19.  Then  the 
transition  from  a  pulsed  to  a  steady  state  continuous  uniform  flow  was  made 
In  section  3.2  and  the  rates  of  molecular  incidence  were  determined  from 
equations  3.2.1  and  3.2.2.  Particular  cases  appear  In  section  3.2.2.  This 
completes  Phase  I  of  the  analysis  of  the  project. 
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4.0  SINGLE  SURFACE  DISCONTINUITY 

A  surface  discontinuity  is  established  by  either  a  discontinuity  in 
the  molecular  flow  or  discontinuity  in  the  pumping  coefficient  or  both. 
Specifically,  the  function  representing  a  discontinuous  flow  or  a  discontinu¬ 
ous  pumping  coefficient  will  be  in  the  form  of  a  step  function.  Further,  the 
discontinuous  line  or  lines  on  a  surface  will  be  a  circle  or  parallel  circles 
on  the  surface.  One  line  discontinuity  will  be  referred  to  as  "single  sur¬ 
face  discontinuity"  and  many  line  discontinuity  as  "multiple  surface  discon¬ 
tinuity".  Except  for  these  discontinuities,  the  molecular  flow  as  well  as 
the  pumping  coefficient  will  be  assumed  uniformly  distributed  on  each  surface. 
The  spherical  configuration  of  the  chamber-vehicle  system  and  the  types  of 
discontinuities  considered  give  rise  to  an  axiosymmetric  case  which  simpli¬ 
fies  somewhat  the  mathematical  analysis  of  the  problem.  As  a  result  of  this 
axiosymmetry,  it  is  intuitively  true  that  the  distribution  of  molecular  hits 
on  a  surface  x  depends  only  on  the  angle  0  between  the  axis  of  symmetry  and 
the  radius  vector  from  the  co-center  of  the  system  to  any  point  on  the  sur¬ 
face  x. 

The  following  analysis  deals  with  a  single  surface  discontinuity  on 
the  chamber  wall  separating  two  spherical  caps.  The  molecular  flow  and  the 
pumping  coefficient  in  each  cap  are  uniformly  distributed  but  differ  from  one 
cap  to  the  other . 

4.1  ANALYSIS  APPROACH  TO  THE  PROBABILITY  EQUATIONS  OF  MOLECULAR  TRANSFER 

The  mathematical  analysis  of  the  single  surface  discontinuity  is  divided 
into  two  parts.  Fart  I  deals  with  setting  up  the  various  probability  equa¬ 
tions  of  molecular  migration  from  one  surface  to  specified  regions  on  the 
same  surface  or  on  another  surface.  These  regions  are  spherical  caps  bounded 
below  by  a  latitude  9^.  Then  the  probability  that  a  molecule  of  a  surface 
will  hit  a  spherical  zone  on  either  surface  will  be  the  absolute  value  of  the 
difference  of  two  probabilities  associating  the  point  p  with  two  spherical 
caps.  Accordingly,  there  are  three  cases  discussed. 

The  first  case  is  when  the  point  p  is  considered  on  the  chamber  wall 
and  the  specified  region  on  the  vehicle. 

The  second  case  is  when  both  the  point  p  and  the  specified  region  are 
on  the  chamber  wall. 

The  third  case  is  when  the  point  p  is  considered  on  the  vehicle  and 
the  specified  region  on  the  chamber  wall.  Thus,  Part  1  will  supply  us  all 
the  information  needed  to  proceed  to  Part  II. 

Part  II  deals  with  computer  solutions  of  the  problem.  Both  surfaces 
are  divided  into  a  sufficiently  large  number  of  spherical  zones  and  the 
molecular  flux  density  on  each  spherical  zone  on  both  surfaces  is  determined 
at  any  time  interval.  The  purpose  of  these  solutions  is  to  provide  us  with 
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sufficient  information  so  that  we  can  determine  the  molecular  flux  distribu¬ 
tion  and  the  incidence  rate  at  various  parts  of  the  vehicle  and  chamber  wall 
as  a  result  of  two  different  steady  state  uniform  flows,  one  introduced  at 
one  part  of  the  discontinuous  surface  and  the  other  on  the  other  part  of  the 
discontinuous  surface. 


4.1.1  From  tho  Chombor  Wall  to  tho  VohicU 


The  transformation  of  coordinates  from  one  orthogonal  system  to  another 
can  be  accomplished  by  means  of  the  following  table 


where 
a,  "E, 
while 


i  j  K 


1  5x 
ha  5a 

1  5v 
ha  5a 

1  52 

ha  5a 

1  5x 

1  Sy 

1  hz 

hp 

hp 

5p 

hp 

sp 

1  5x 

I  5v 

_1 _ ^ 

hy 

dy 

h7 

S7 

h7 

57 

(A.l) 


JL,  j,  K  are  unit  vectors  in  the  direction  of  x,  y,  z  in  one  system  and 
c  are  unit  vectors  in  the  directions  of  01,  p  and  y  in  the  other  system, 
ha,  hp,  hy  are  defined  as: 


hoe 


Ss 


hy 


hi 

^7 


Now,  let  P  be  a  point  on  the  chamber  wall  with  position  <I>p  and  6^  and  let 
be  taken  in  the  directions  of  increasing  G  and  r^at  P.  ^Then  y 
coincides  with  the  normal  to  the  surface  at  P  drawn  outwards,  and  the  plane 
(a,  p)  becomes  the  tangent  plane  at  P.  The  element  of  length  on  a  sphere  of 
radius  r  is  given  by 


2  2  2  2  2  2 
ds^  »  r  cos  0d^)  +r  d0 

2  2  2 

The  element  of  length  at  P  is  given  by  ds  «  dr  +ds-  where  ds-  is  an  arc 
element  on  the  sphere  passing  through  P. 

Hence 


ds‘ 


.  2^  2  2^  ,^2^  2  .^2 
dr  +r  cos  6d^>  +r  d0 
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and  \CC  ■  rco80,  hp  ■  r,  h7  •  1 

follow  immediately. 

The  equation  of  a  sphere  is 


X  -  r 


|co8^cos0i+sin®cos0j+sin0^ 


Then  on  the  sphere  at  F 

Sx/Soi  *  -r  sin®co80  Sy/Sx  ■  r  co8®cos0 

3x/Bp  =  -r  cos$sin0  5y/Sp  *  «r  sinC>sin0 

Sx/Sy  =s  cos$cos6  hy/^y  =  sin$cos6 

5z/Sa  =  0,  Sz/5p  =  r  COS0,  Sz/S  =  sin0 

and  the  table  (4.1)  becomes 

i  j  K 

a  -sin(I»  cos$  0 

"B  -cos$sin6  -8in08in6  cos0 

c  cos^cosd  sin<l^cosd  sin6 


(4.2) 


As  a  result  of  axiosymmetry,  the  probability  that  a  molecule  will  hit  the 
spherical  zone  (d . ,  d.)  is  independent  of  Thus  we  may  set  ^  »  0  in  the 
table  (4. 2) to  obtain  ^ 


(4.3) 


In  order  to  preserve  the  notation  used  previously  in  connection  with  this 
project^  a  will  denote  the  angle  between  the  radius  vector  from  P  and  *-c  and 
P  will  be  the  angle  between  a  and  the  projection  of  the  radius  vector  from  P 
on  the  plane  (a,^). 


i 

j 

K 

0 

1 

0 

-sin0 

P 

0 

cosS 

] 

cosd 

p 

0 

sind 

1 

Then  (see  Figure  1)  by  means  of  table  (4.3)  we  obtain 

X  ■  r  sin^cosd  a 
o 

+r  [-cos^coad  sind  +8in8  cosd  b 

L  op  ®  _ 

4t  fcos^cosd  cosd  +sind  sinP  1  c 

L  o  P  o  pj 


(4.4) 
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■  Rcos/  jcospsin7a+stn3sin7b-cos7cJ  (A.  5) 

Now  let  the  heads  of  vectors  x  and  y  coincide.  They  y  »  PO  +  X  which  yields 

(4.6) 


y  B  r  sin$cos6  a 
o 


-c 


+r  -cos<I>cos0  sin0  +sin0  cos6  I  b 
op  op 


+  <r fcos<I>cos0  cos6  +sin6  sin6 
\  op  op 


The  intersection  of  C  and  0^  can  be  found  if  we  set  ~  V 


Then  we  obtain 


a)  sin$cos0^  =  cos7Cosp 

b)  -cos<l>cos0  sin0  +sin0  cos0  =  cosysinB 

o  p  o  p 

c)  cos$cos0  COS0  +sin0  sin0  =  sin7 

o  p  o  p 


(4.7) 


from  4. 7c  we  obtain 


COS0 


— — —  -  tan©  tan0 
COS0  cosd  o  p 


(4.8) 


substituting  4.8  into  4.7b  and  solving  for  sin^  we  obtain 


sinp  = 


sin0  -sinysin© 

O  ' _ £ 

COS7COS0 


(4.9) 


If  we  take  the  dot  product  of  the  vectors  y  (eq.  4.6)  and  PO  we  obtain 

2 


2 

cos  a 


^R*»r  |co8(i)cos0^cos0^+sin0^sin0^^ 


|y|' 


But  from  4.6  we  have  also 


|y|' 


H 


1+sin  7-2sin7(cos$cos0^cos0p+sin0^sin0p) 


Hence 


2  2 1 

o  sin  7  l-(cos^cos0  cos6  4-sin0  sin0  ) 

.  _  o  p  o  p 

sin  a - 


1+sin  7-2sin7(cos(I>cos0p+sin0^sin0p) 


(4.10) 
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from  equation  4.6  we  obtain 


r  8ln<I>co8d 


|y|co8^8ina 


Hence 


-cos^cosd  sind  +8in6  cos6 
op  op 


|y| 8ln38iTa 


tan^ 


-co8<I>cosd  sin6  +sin0  cosd 
o  p  o  p 

sln^coed 


differentiating  4.11  with  respect  to  ^  we  obtain 


d^  * 


s 


cosd  Sind  cosd  -cos$sind  cosd 


l-(cosC)cosd  cosd  +sind  sind  ) 
o  p  o  p 


Summarizing  the  above  results: 


The  points  of  intersection  of  the  curves  and  d  are  given  by: 

®p 


cos$  * 


siny-sind  sind 

-  _ Q _ P 

cosd  cosd 
o  p 


in  terms  of  $ 


sinp 


8 ind s iny s ind ^ 
cos/cosd 


in  terms  of  p 


On  d 


sin  a  - 


sin  y  l-(cos$cosd  cosd  +sind  sind  ) 


1+sin  7-2sin7(cos®co8d^co8dp+8ind^slndp) 


and 

cosd 


dp 


sind  cosd  ••cosOsind  cosd 
_  p  o  o  p. 


l-(cos$cosd  cosd  +sind  sind  ) 
'  o  p  o  p^ 


(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 
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Let  max  $  niin  denote  the  maximum  and  minimum  latitudes  of  curve  which 
p  P  p 

occur  at  $  ■  0  or  ®  »  jc. 


When 

-f<0p<-7 

from 

Eq.  4.8  we  obtain 

at 

o 

1 

max  C^ 

P 

at 

min  C^ 

P 

■-<!  +  !■  +  «p) 

When 

-7<®p<7 

from 

Eq,  4.8  we  obtain 

at 

4>  =  0 

max  C^ 

P 

■ !  ^  ®p  ■ 

at 

$  =  0 

min  C^ 

P 

-  -f  +  +  p 

When 

again  from  4.8  we  obtain 

at 

$  “  K 

max  C^ 

P 

=  f  +  r  ■  Sp 

at 

a.  -  0 

min  C^ 

P 

-  -  f  +  r  +  «p 

Due  to  the  geometrical  symmetry  it  will  suffice  to  consider  the  case  ^^>0. 

We  distinguish  two  cases 

Case  I  ^6  >  max 

2^  o-  6 

P 

Case  II  max  >6  ^  min 

6  o 
P 


Case  I  - 

e  >7.  f  >  max  Cg 
^  P 

(see  Figure  2)  6  does  not  intersect  C^  •  Let  be  the  surface  region  on 
the  vehicle  bounSed  below  by  6  •  Let  P(§  ,6  )  denote  the  probability  that  a 
molecule  from  the  point  p  will^hit  the  region 
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Then  with  the  aid  of  4.15  and  4.16  we  obtain 


P(0  »6-)  ■  “//  sinacosadad3 


P  o 


slnd  cosd 
P  o 


jln0  siny  2*  "  sing  cos©  ^  ^ 

ry  f  ^  rt 


4k 


O  1J-  -I  ^ 

cos®  +tan0  tan0  .f— S,  ,7 


o  P  2sin7cos6^cos0p 


P(0  >0  ) 

p*  o 


Or 

sine  sin/ 


(1+sin  y)sine^-2sin7sine^  2jc 


djccose  cose 
o  p 


I 


d$ 


(4.17) 


o  cos®+tane  tane  -  l+8in27^ 


o  p  2sin/co8e^co8ep 


In  order  to  evaluate  the  integral  in  4.17  we  use  the  residue  theorem  by  setting 
2co8®  -  Z+Z’^ 


Let 


A  “  tane  tane 


l+sin^y 


o  p  Zsin/cose^cosep 


Then 


L 


2k 


d® 


o  cos®+A 


lx _ ^ 

^  Z^+27^+1 


(4.18) 


The  roots  of 


Z  +2>^Z+1  -  0  are 


(4.19) 


-?d:(7^^-l)^^^ 
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Consider  the  equation 


we  obtain 


f(x)  -  x^+2?sfl 


f(l)  -  l+2>^+l  -  2(l+>s) 
f(-l)=  1-27V+1  .  2(1-7^) 

Asinycos  — — ^  +  (l-siny)^ 

1-A  *  - ? _ _ 

2sin7Cose ^cose^ 

-4sin7sin  P  -(l-siny)^ 

l+T,  - - — - 2 - - 

2sin7cose  cos© 
o  p 


(4.20) 


(4.21) 


2  9 

Hence  1-^  <0  or  A  -1>0 

conclude  roots  and  on  account  of  4.20  and  4.21  we 


-  K-V  (X^- 1)  ^/^<1<-?,+(7^2_  ^  j  1/2 


Hence  the  residue  of 


Z  +2?v2+l 


within  /Z/«l  is 


2  (7v^- 1)^/2 


it  follows  2jt  , 

f  —  -2rt 

cos^A  ,^2  ,,1/2 


(A  -1) 


Therefore 


where 


P(e  ,0  ) 

p*  o' 


sine^siny  (1+sin  7)sine  -Zsinysin^ 


4cose  cose 
o  p' 


(4.22) 


sin0  siny  _ 

o  ^  ^ 

2  "  2 


Kl+sin  7)8ine  -Zsinysine  1  siny 


I —  o  Pi  _ 

£2sin7Cos(e^+ep)+l+sin^^  -28in7co8(e^-ep)+l+sin^7 
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we  can  easily  verify  that  when  ®  ®o  with  (i,e,6^  =  y 

2  °  P  ° 

that  P(e  ,6  )  ■  sin  y  as  it  is  expected 
p  o 


When  e  -  e 
o  p 


sin6  sin/  (1- sin/) sin/s in© 


2  2  1/2 
2(l+sin  /+2sin/cos  6^)  ' 


■  2 


Case  II 

6  always  intersects 
o  6 


This  occurs  when  max  >Q  >min  whether  0  >/  or  0  </  see  Figures  3. 

00  0  p'  p'  ° 

P  P 

The  intersection  of  0  and  is  given  by  cos$  =  - 2“  "  tan0  tan0 

00  COS0  COS0  o  p 

p  o  p  ^ 


4.23 


(see  Eq,  4.8).  Let  $  be  the  first  intersection,  i.e.  0<<I)<it.  Then  the  other 
intersection  is  °  Let  be  the  surface  area  on  the  vehicle  bounded  by  0^ 

and  and  for  wh?ch  0>0  .  p 
0  o 

P 


Then 


p(e  ,e  )  =  -  //  siiciicosadadp  = 


sin0  COS0 
-2 _ 2 


sinews in/  o  (cos$-3in^^-03ep  ) 

°  1+sin^y 

cos<I4-tan0  tan0  -  ^  .  * 


sin^y 


o  p  2sin/cos0^cos0p 


where  is  given  by 


sinp 


sin0  -sin/sin0 

_ _ -O _ 


O  COS/COS0 


(Eq.4.9) 
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or 


sin©  8in7 

■  — - 


(4.24) 


(l+sin  7)sin0^-2sin7sin0^ 

4jtcos0^cos0 
o  p 


0) 

;  ° 


d(ii 


cos4>+tane  tan©  -  ^  - — 

o  p  Zsin^cosd  cos6 


+  f  -  3,) 


To  evaluate  4.24  use  the  following  substitution 


X  =  tan  2 


d<t 


2dx 

1+x^ 


cos$  =  2cos^  4  -1  =  ^ 

14x 


(4.25) 


Then 


d$ 


josflH-A  >d-l 


h-  ! 


dx 


2  X+1 


'VI 


where 


Now 


X  =  tan©  tane  -  ^  7 - - 

o  p  2sin7Cos6^cos6p 


^  2  >0  since  A^-1>0 

(A-1)^ 


Solving  4.25  for  X  and  setting  $  »  we  obtain 


co8(0^*0p)"sln7 


CO  s  (6  ^+6  p  )+s  in7 


1/2 


where  we  have  taken  the  positive  root  since 


(KO  <jt- 
0 


►  tan  -  X^O 


22 


AEDC.TDR-63«88 


Then  4.24  becomes 

Sind  8ln7  .  2  . 

p(e  ,e  )  -  — :r -  ®  °  ?  -p  )-G(e  ,e  )tan" 

'  p*  o  2jt  o  jc  '  2  ^  p*  o 


‘i(M) 


1/2 


(4.26) 


where  G  has  been  defined  in  Eq.  4.22 


and 


-  -  28in7COs(6  +0  )+l+8in  7 
2^  .  _  _  -O  P _ ' 

>+•1  -2sin7Cos(e  -6  )+l+8in'^7 
o  p 


The  formulas  4.22  and  4.26  should  have  the  same  value  in  their  common  region 
of  domain. 


Set 


e  ■  Max 
o  6 


in  Eq.  4.26 


and 

Then 

Hence 


A  -  i  Q  -  21 
®o  2  *  ^o  2 


->00 


tan 


-1 


*i(M) 


1/2 


2cl 

« 

2 


^  0 


and  4.27  reduces  to  4.22. 


Also  when  6  =  min.  C,  we  expect  P(0  ,6  )  •  sin  7 

o  G  '  p'  o  ' 


Indeed,  0=-“+0  +  7 

o  2  p  ' 


Xi  =  o- 


•  tan 


-1 


^  «  0,  p 

o  o 


% 

’2 


-.1/2 


2ri 


7&l 


»  8in"7 


Consider  the  special  case  where  6  ■  0  *  0 

op 


Then 

Then 


sinPo  -  0  Po  -  0 


(what  was  expected) 
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4J.2  From  th«  VthicU  to  tho  Chombtr  Wall 

Consider  Che  point  P  on  the  vehicle  and  the  curve  6  on  the  chamber  wall 
(Figure  4,3).  Again  we  need  consider  only  ^^^0  with  two  Bases, 


Case  I 


■^>  max  Cg 


Case  II  .Max  C^  >6  >  min  C- 
o  "  cr“  t7 


(Figure  5) 
(Figure  4) 


Let  P  be  a  point  on  the  vehicle  with  position  (0,  6^)  and  let  A  be  a  point  on 
6^  such  that  01^  (see  Figure  4). 

If  we  set  up  an  orthogonal  system  at  P  as  before  and  make  use  of  Table 
(1,3)  we  obtain 


X  «  RsinOcosd  a 
o 


+  R 

+  R 


-cos^cosd  sin0  +sin0  cos0 
op  op 


cos^cosd  cos6  +sind  sin 
op  op 


] 


■] 


(4.27) 


y  ■  Rsin^cosd  a 
o 


-c 


+  R  -cos^>cos6  sin0  +sin0  cos6 
o  p  o  p 


(4.28) 


+  Rjfcos^tcos©  cos6  +sin0 

IL  °  p 


sin6  l-A 
O  pj  ) 


The  intersection  of  the  curves  C^  and  6  (if  they  intersect)  is  easily  found 
to  be  given  by  p  ° 

cos4>cos8  cosd  +sin8  sin0  *  1 
op  op 


Hence 


as  before. 


cos^>  “  .SilS - -  .  tan0  tan6 

o  cosd  cos8  o  p 


(4.29) 
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Also 


"  2  I  2  ” 

|y|  ■  R  1+sin  /-Zsin/Ccos^cos^^cos^^+sin^^sinSp) 

2 

l-(cos<)cos0  CO80  +sln6  sin0  ) 

Sin^a - 5 - 2 - E - 2 - 2 - 

1+sin  7-2? In7(co8®co86^cos6p+8ine^8inep) 


(4.30) 

(4.31) 


sin3 

o 


sinO  -sinysinO 
o  ^ _ P 

cosG  cosy 
P 


(4.32) 


dp 


cosd  slnd  cos6  -cos<&sin6  cosO 
oL  P  o  _ o  p. 


d$ 


l-(cos4>cos6  cos6  4-sin0  sinfi  ) 
o  p  o  p 


(4.33) 


In  an  analogous  manner  as  before  (sec.  4.1.1)  the  following  relations  can  be 
derived 

When 

at  $  =  jt  Max  Cq  =  t  -  0  +7 

ep  2  p 

=  0  Min  =  -~  +  e  +7 

C7  2  p 

p 

When  ”7<^p<7 

at  (I)  =  0  Max  =  “  +  6  -7 

v  4  p 

P 

$  =  0  Min  =  -~  +  0  +7 

(7  4  p 

P 


When 

at 


- 

$  =  0  Max  Cg  *  2  ^ 

P  ^ 

®  -  «  Min  Cg  »  -(|  +  e  +  7) 

T\  ^ 
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Case  I  r>e  >  Max  (Figure  5) 

—  cr~  u 

P 


which  d>6  a 
—  o 


It  occurs  when  6  >7.  Let  be  the  surface  area  on  the  chamber  wall  for 
^  P  “  ^ 


Then 


p(e  ,e  )  -  -  //  siiacosodadp 


p  o 


(4.34) 


slnd  cosO 

(cos$  -  -st^PcosAp  '>  ^ 


°  l+sin^y 

cos^tanfi  _  tani0  -  ^  ^  ^ 


o  p  2sin7cos0^cos0p 


sin0^  (1+sin  7)sln0^"28ln7sln0^  2jt 


do 


2sin7 


dicsln  7Cos6  cosd 
o  p 


cosO+tane  tan©  -  1+  sin  7 - 

o  p  2sin7CO86^co80p 


Sind 


o  icG 


28ln^ 


Consider  the  special  case  where  0  “  2*  ^o  expect  P(6  ,  0^)  *  !•  A 

substitution  easily  verifies  that^this  Is^indeed  the  case.  ^  ^ 


Case  II  Max  C^  >0  >Min  C^  (Figure  4) 
y  O  a 

P  P 


As  before  we  find 


sin0  . 

P-O^  -  Po> 


(4.35) 


.  -1 


P  o 

,  2 
sin  7 


tan 


*i(^) 


1/2 


where  G  is  defined  In  the  same  way  as  for  equation  4.26. 
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of  domain 


Foraulas  4.34  and  4.35  should  slve  the  same  value  in  their  comnon  region 
lain.  Indeed  Ife-HaxC.  -f-S  +  7  6  >y 


*0  "  ^  "  2* 


P 

1/2 


n 

2 


4.35  reduces  to  4.34* 


Also  when  6  ■  Min  we  expect  P(e  ,0  )  -  1 

O  C7  P  O 

P 


Set  •  -f  ®p  +  7»  7<0<  f  or  -7<B^<7 

in 


»o  ■  ^  -  !• 


and  the  result  follows  immediately. 


‘i(M) 


4.1  *3  From  tho  Chambor  Wall  to  tho  Chombor  Wall 

Consider  both  the  point  F  and  the  ctirve  6  on  the  chamber  wall.  Thus 
(see  Figure  6) 

y  *  -Rc  +  X 


Rsin^cosd  a 
o 


(4.36) 


+R  I-costcosS  sin©  +8lne  cosS  b 
o  p  o  p 


L- 


-fR  cos$cos0  cos6  +sin0  sin0  -l|  c 
o  p  o  p 


-1 

3 


y^«  2Rcos7  l^ospslny  a  +  sinpsin7  b-cos7  cj 
For  the  intersection  of  6  and  (if  they  do  intersect)  we  set 

O  V 

-  -  P 

y  ■  y.  Then 


(4.37) 


a)  slittcosd^  ■  2co87C083sln7 
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b)  -cos<I>cos6^sin6p+sin6^co80p  »  2co878lnP8in7 

2 

c)  cos<I>cos0  COS0  +81110  8in0  -1  •-2co8  y 

o  p  o  p  ' 


from  4.38c)  we  get 


co8(I>  =  '  cos27 —  _  tan0  tan0 

COS0  CO80  o  p 

op 


Substituting  4.39  into  4.38b  we  get 

sind 

sinp  =  cotZytanep  + 


Similarly  we  find 


2  1 

sin  a  =  “r  (l+cos(I>cos0  cos0  +sin0  sin0  ) 
2  op  op 


COS0 


d3 


0_L 


sin0  COS0  -cos<I>sin0  cos0 


l-(cos4>cos0  COS0  +sin0  sin0  y 
o  p  o  p 


f  -  27>0 

f-27>ep>27-f 
Max  Cg=  27-ep 

Min  Cg=  -(6+27) 
P  ^ 


f  -27<0 


27-^  -27 

Max  C.  »  X-  (27-0  ) 

0  '  '  p' 

P  ^ 


Min  C 


0 


-X  +  ly^ 


2y-fe>  -  f 


Max  C 


0 


it-(27-ep) 


f  -27  >e  >  -  f 

2  P  2 

Max  -  x-(27-0  ) 

0  P 

P  ^ 


Min  =  -(0  +  27) 

0  p  ' 


Min  *  -(27  +  0  ) 

0  '  p 

P 


(4.38) 


(4.39) 


(4.40) 


(4.41) 

(4.42) 
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““  Ce  -  2r-«p 

Min  C-  -  -t+2y+e 
%  ^ 

Again  we  distinguish  two  cases.  Case 

Case  II  when  6  does  intersect 

0  u 


f^>27  -  f 

Max  -  27-e^ 

Min  C.  -  -*+2yfe 
®P  P 

I  when  6  does  not  intersect  and 
o  y 

p 


Case  I  Max 


P(e  ,0^)  ~  ^  If  sinOcosadadp 


where  is  defined  as  previously  but  on  the  chamber  wall.  It  is  easy  to  see 
then  that 


■  h  4 


2ic 


cosd 

sind  cosd  -cos^sind  cosd 

d« 

0 

P  0 _ 0  P- 

l+(cos^cos6  COS0  cosO  +sini9  sin9  ) 
o  o  p  o  p 


j  sin9  COS0 

8in9  2jc  8in0  cos0 

_ o  r  \ _ o  & 


:) 


cos^  tan0  tani9  +  '  ^  ^ 

o  p  cosd  coaO 

•  A  r 


Sind  Sind  sind 
_ o  , _ o 


2x 


^xcoae^coaSp  eo8*+tane_tan0_  + 


O  ~  p  CO80  CO80 

op 


Let  "K  ■  Can0  Caii9  + 


O  p  CO80  CO80 

op 


,  1/2  8ine  +8in0 

a  -1)  -  >0 


CO80  eos0 
o  p 


a 
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As  before 


f(l)  -  2(1+A) 


f(-l).  2(1-7.) 

cos (6  -6  )+l 


f(Z)  -  Z  +2>>Z+1 


(!+>.)« 


o  P 


2co8^  ^2;^ 


cos6  cosd 
o  p 


cosd  cosd 
o  p 


>0 


cos(0  +0  )-l 

1-^ - ■  2  P-  <0 

COS0  CO80 
o  p 

2  1/2  2  1/2 
Hence  .^(X  -1)  <-l<-X+(X^-l)  <1 


and  the  residue  of 


Z  +2XZ+1 


within  IZl  -  1 


is 


2(7.^-!)^/^ 


It  follows  /  - -  7  # 

■'o  cosW-7.  i  „2 


Z  +27.Z+1 


2ii 


therefore 


,  sind 


(4.43) 


Case  II  Max  C.  >6  >Min 

0  0  6 
P  P 

sin9 


"VV 


*i(M) 


1/2 


(4.44) 


30 


AEDC.TDR.63-88 


Where 


®  *  cos 

o 


-1 


*^008^7+81110  8in0 
^  o  p 


P  »  sin 
o 


-1 


CO80  COS0 
O  p 


8in0 


cot27tanep+ 


CO8(0^-0  )  +  C0827 


A-1 


A+1 


cos(0  +0  )  -  C0827 
'op  ' 


1/2  l-cos(0  +0  ) 
_ '  o  p 

*  sin0  +siii0 
o  p 


1/2 


Notice  that  equation  4.34  reduces  to  equation  4.53  when  0  ■  Max  where 

mM _  ^  _ rt-..  ^ _ _  _  ® 


Max  “27-  0  .  Since -then, 

®p  p 

*o  "  ^  “  2* 


1/2 


It 

2 


When  0  *  0  *  0,  then  $  «  «-27,  p  *  0 

p  o  *  o  '  o 


tan 


-1 


*i(M) 


1/2 


P  o 

4.2  RECAPITULATION 


and  equation  4.54  givea  us 
for  the  half  illuminated  sphere. 


A.  From  the  Chamber  Wall  to  the  Vehicle 

Extending  the  previous  results  to  the  whole  chamber  wall  we  define 
as  the  surface  area  on  the  vehicle  which  is  accessible  to  the  point  P  ^ 

and  for  which 

0>0 
“  o 

and  P  (0  ,0^)  as  the  probability  that  a  molecule  leaving  the  point  P  will  hit 
the  region  . 

P 


(a)  0 

o 
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P(e  ,e  ) 

'  p’  o 


sine^sin/ 

2  "2 


(b)  Min  C-  <0  <  Max  C 


sln0  sin/ 


where 


p(e  ,e  )  -  , 

'  p’  o  2« 


®  + 

o  « 


,  2 
sin  y 


(!  -o) 


-G  tan 


-1 


1/2 


-2sin7C08(6^-6p)+l+sin^7 


cos(0  -6  )-sin7 

V  ^  op'  ' 

1/2 

1  cos(d  +6  )+sin7 

|(l+sln^7)sln0^-2sin7sln9^ 

sln7 

jt  ‘^^sin7COs(e^+6p)+l+8in^ 

”  2 
-28in7CO8(0^-0p)+l+8ln  7 

and  are  given  by 

: 

r 

“1 

1/2 


$  =  cos 

o 


- ^ - T  -tan0  tan9 

COS6  COS0  o  p 

op  ^ 


,  sin0  -sln7sin0 
p  .  sin'^  - 2_-- - E 

O  COS7COS6 


with  the  restrictions  that 


0<a>  <ic 

o 


■!«oi 
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(4.45) 


(4.46) 


(4.47) 


(4.48) 


(4.49) 


(4.50) 


(4.51) 
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B.  From  the  Vehicle  to  the  Chamber  Wall 


Now  e  and  9  switch  surfaces.  is  defined  as  previously.  is 

°  P 

the  intersection  of  the  tangent  plane  at  P  (of  the  vehicle)  with  the  chamber 


wall. 


(a)  G  completely  illuminated  from  P 

^0 

sin0  - 

P(0  )  -  TT^  " 

0’  o  2sin7  28inS 

(b)  Min  C_  <e  <  Max  C. 

t?  cr~  V 


(4.52) 


slnd 


p*  o  2itsin7 


1 

/fl  _  \ 

G  -1 

(1>  +  - 
0  It 

[T^o] 

-  — j  tan 
sin  y 

*1  \7&1 ) 

1/2 


(4.53) 


where 


are  given  by  4.47,  4.48,  4.49,  4.50  and  4.51  respectively. 


C.  From  the  Chamber  Wall  to  the  Chamber  Wall 


Both  9^  and  are  on  the  chamber  wall.  (now  on  the  chamber  wall) 


0  p 

is  defined  as  previously 


(a)  9  >  Max  C^ 

'  o  9 


sin9  . 

- 2  +1 

2  2 


(4.54) 


(b)  Min  Cg  <e^<  Max  Cg 
P  P 

sin6 


where 


+  7  (r^o) 

1/2 


2  .  1  -1 
—  tan 

1C 


‘i(M) 


1/2 


(4.55) 


are  given  immediately  after  equation  4.44. 
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5.0  COMPUTER  LOGIC 


The  problem  of  determining  analytically  the  molecular  flux  distribution 
in  the  cases  where  discontinuous  out gassing  characteristics  have  been  assumed 
is  exceedingly  difficult.  Consequently  we  have  used  a  digital  computer  to  ob¬ 
tain  solutions  to  the  problem.  The  types  of  surface  discontinuities  we  deal 
with  have  been  explained  previously  (see  4.0  Single  Surface  Discontinuity). 

The  procedure  which  is  followed  for  a  computer  solution  consists  of  parti¬ 
tioning  the  surfaces  into  a  number  of  spherical  zones,  each  spherical  zone 
having  uniform  outgassing  and  pumping  characteristics.  The  probability  of 
molecular  transfer  from  each  point  of  a  spherical  zone  to  another  spherical 
zone  will  be  considered  constant,  thus  giving  rise  to  a  step  function  along 
the  meridian  of  a  surface  for  the  probability  of  transfer  from  a  point,  to  a 
spherical  zone.  The  molecular  density  of  the  molecules  which  leave  a  spherical 
zone  will  be  taken  as  uniformly  distributed  on  this  spherical  zone  so  that  the 
molecular  density  along  a  meridian  will  assume  the  form  of  a  step  function 
also.  Thus  the  exact  solution  of  the  molecular  flux  distribution  can  be  better 
approximated  by  increasing  the  number  of  the  spherical  zones  of  each  surface. 

5.1  PARTITION 

Both  surfaces  (chamber  wall  and  vehicle)  are  partitioned  into  2N 
spherical  zones  of  width  t£  where  N  is  an  integer.  In  numbering  the  sequence 
of  these  sperhical  zones  we  set 

\  (r-N)Z^] 

where  has  been  defined  previously. 

*1  ■  [-f  •  -  f  *  '^] 

Xj  -  [(i-H)ie,  (2-ii)ie] 

W  -['>•  *] 
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The  choice  of  N  is  such  that 

a)  N  is  reasonable  large. 

b)  The  discontinuous  line  or  (lines)  are  boundaries  of  spherical  zones. 

c)  G  »  y  is  a  boundary  of  a  spherical  zone* 

Obviously,  the  probability  from  a  point  of  a  latitude  of  a 
surface  X  to  a  spherical  zone  is 


(5.1.1) 


In  terms  of  the  partition  then  Eq.  (5.1.1)  will  denote  the  probability 
from  a  spherical  zone  X  to  a  spherical  zone  Yj  producing  therefore  a  step 
function  for  the  probability  distribution  of  a^ surface  with  respect  to  a 
spherical  zone*  Clearly  then,  a  sufficient  condition  that 

P  »  p 

Vr 


is  that 

K  =  2N-I+1,  r  =  2N-J+1 

Hence 

Pv  V  •  Pv(2N-I+l)  Y(2N-J+1)  (5.1.2) 

Hence,  in  view  of  Eq.  (3. 3. 2)  it  is  sufficient  to  calculate  P  for  J>  lH-1 


5.2  PROBABILITY  EQUATIONS 


The  probability  equations  which  will  be  used  in  the  computer  solutions 
have  been  derived  previously  and  appear  in  section  (4.2).  In  regrouping 
these  equations  for  computer  use  we  have 


given  by  (4.45) 
given  by  (4.46) 


from  chamber 
wall  to 
vehicle 


ww 


i-sine. 

2 

l-sin0j 

~2 


sin^7 


from  chamber 
wall  to 
chamber  wall 


Pcc(^I>^j)  “  given  by  (4.55) J 


(5.2.1) 

(5.2.2) 

(5.2.3) 

(5.2.4) 

(5.2.5) 
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Equations  (5. 2.1)  and  (5.2.2)  pertain  to  the  case  "From  Chamber  Wall  to 
Vehicle"  while  Eqs.  (5. 2. 3)- (5«2, 5)  pertain  to  the  case  "From  Chamber  Wall 
to  Chamber  Wall".  Here  we  have  omitted  the  case  "From  Vehicle  to  Chamber 
Wall"  in  view  of  the  relation 

^VC^^I*^J^  ”  ^  (5.2.6) 

Thus  there  are  actually  two  probability  matrices  which  have  to  be  determined, 
one  for  the  chamber- vehicle  case  and  the  other  for  the  chamber-chamber  case. 

5.3  PROBABILITY  FLOW  CHARTS 

In  all  cases  we  determine  the  probabilities  (P„Y(e-,8-)  for  J>  N+1.  In 
differencing  the  results  by  means  of  Eq.  (5.1.1)  ana^excenaing  the  domain  of 
J  to  J<IW-1  by  means  of  Eq.  (5*1.2)  we  obtain  the  probability  matrices 

•”1  fee- 

5.3.1  From  Chamber  Wall  to  Vehicle 

Since  it  is  sufficient  to  start  with  view  of  the  relation  (5.1.2) 

we  may  begin  with  0^  y~it/2.  For  each  ^jWe  have 

Lim  C^  =  0_+Jt/2-7  if  0_-7<O 

=  -0j+jt/2+7  if  6j-7>0 

As  for  the  L^m  C  ,  we  are  only  interested  in  the  value  of  Lim  C^  when 
Lim  C.  >0.  ^I  This  occurs  when  e_>rt/2-7  yielding  ^I 

Lim  C.  =  e^+y-n/2 
For  0j-7<O 

the  domain  of  6  for  a  non  zero  probability  is  0<6  ^<B  -hc/Z-y  Lim  C„  <0 
J  J  1 

Lim  C^  >0 

and  i-s  determined  from  Eq.  (5.2.2) 

For  ej.-7>0 
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the  domains  of  6j  are 

a)  q^j<-ej-hi/2+7 

or  0^47-it/2^j<-0j+n/2+7 

b)  -0j+«/2+7^^N 


Lim 

Llm 


<0 

I 

>0 

I 


and  P^^(0^,0.)  is  determined  by  Eq.  (5.2.2)  for  case  a)  and  by  (5.2.1)  for  case 
b).  In^the  flow  chart  we  have  been  using  spherical  zone  sequential  numbers. 

The  corresponding  angles  to  these  numbers  are: 

N-N2+I - ^0  -  n/2-7f0j  -  Lim  0j-7<O 

IH-N^+l - ^0  -  7 

2N-N2+I - ^n/2-y 

N^-N+I - ^0  »  7-jc/2+0_  -  Lim  C. 

2  I  0j 

3IWJ2+2-I - ^0  -  n/2+7-0j  -  Lim  0j-7>O 


5.3.2  From  Chomb«r  Wall  to  Chamber  Wall 

We  distinguish  two  cases  I  and  II 
Case  I  n/4-7>0 

We  start  with  0^  “  0  and  Jt/2 

We  have 

Lim  «  0_+n-27  6-<2y~K/2 

01  I  I 

-  27-0^  0^27-k/2 

As  for  Lim  C^-  we  notice  that  Lim  ■-(0-+27)>O  when  -tt/2<0  <-7  and  for  0^-7, 
Lim  <0  in  all  cases.  I 

^I 

a)  For  0j<27-it/2 

the  domains  of  0 ,  are 

if  Lim  C.  >0 
®I 


i)  Q^j<-(ei+27) 

ii)  -(ej+27)^j  e^+K-2y 
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or 


iii)  O^jO^+K-Zy 


if  Lira  C.  <0 


In  case  la,  the  Probabilities  determined  by  Eq.  (5.2.4) 

for  i)  and  by  (5.2.5)  for  ii)  and  by  Eq.(5.2.5)  for  iii) 


b)  For  B^Zy-nll 


the  domains  of 


iv)  0^j<- (6^+27) 

v)  -ie^2y)^^<2y-e^ 


\ 


C  if  Lim  >0 


Vi) 

27-01^  j<n/2 

or 

vii) 

O^j<27-0j 

Lim  C-  <0 

T  —  ' 

viii) 

27-ej.^j<it/2 

or 

J  Lim  Cg  >0 

ix) 

P^j<jt/2 

if  Hm  Cg 

Then  Eq.  (5.2.3)  for  calculating  used  in  the  cases  iv) 

while  Eq.  (5.2,4)  is  used  for  the  cases  v;  and  vii)  and  Eq,  (5.2.5)  for  the 
cases  vi),  viii)  and  ix). 

The  corresponding  angles  to  the  sequential  numbering  which  appears  in 
the  flow  chart  are 


2(1H-1-N.)-I - >e  =  -(6^+27)=  Lim  C. 

i  i  t7- 


2(N-N.)+I - >Q  =  jt-27+e_  =  Lim  C.  ,  Q.<ly-%ll 


2N2+1 — 
N-2N2+1- 


-)-e  =  Zy-ull 


->0  =-27 


Casa  II  n/4-7<0 


In  this  case  we  start  with  6j  =  0  and  0j>27-it 
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-  Sj+Jt-Zy 

=  27-6^ 

As  for  Llm  we  have 
®I 

=  ej.+27-it 


We  have 


Llm  C 


0, 


e^<2y-nl2 

e^ly-n/l 


if  ej+27<K/2 
if  ej.+27>«/2 


c)  e^y-nll 

The  domain  of  9j  is 

x)  O<0,<0_-ht-27  Llm  C-  <0 

J  L 

xi)  0_+27-rt<0 -<0_+j(-27  Llm  >0 
i.  "  J  J,  — —  y 


d)  e^>2y-«/2 

The  domain  of  dj  Is 
xii)  0^j<2y-e^ 

xili)  2y-e^^j<x/2 


Llm  <0 


xlv) 

xv) 


or 

ej+27-n<ej<27-e 

2y-e^^j<x/2 


Llm  C„ 
y 


>0 

I 


Eq.  (5.2.3)  is  used  for 
for  all  the  other  cases  in  c) 


the  cases  xili)  and  xv)  while  (5.2.5)  is  used 
and  d). 
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The  corresponding  angles  to  the  sequential  numbering  which  appears  in 
the  flow  chart  are 


2(N-N2)+I - 

- ^0 

■  6-+K-2y  ■  Lim  C„ 

Sj.  <27-*/2 

2N2+I - 

— >e 

-  27-«/2 

3N-2N2+I - 

— >e 

■  «-27 

2(N2-N)+I - 

-^e 

■  e_+27-«  ■  Lim 

ej+27>K/2 

2(N+N2+1)-I - 

— ve 

“  27-0J  »  Lim  Cg 

e^2y-nl2 

5.4  MOLECULAR  TRANSFER  FLOW  CHARTS 

(single  discontinuity  on  the  chamber  wall) 

The  area  of  a  spherical  zone  of  width  on  a  sphere  of  radius  r  is 

4jtr^sin  ^|cos  (0+^)  (5.4.1) 


If, 

and 


however,  tB  is  small  we  may  set 


sin-2 


2 


the  spherical  zone  area  will  be 
Zjtr^cos  (6+^)Zi0 


(5.4.2) 


Thus,  the  flow  which  will  hit  a  spherical  zone  Y.  due  to  the  flow  which 
leaves  a  spherical  zone  X  is  (Area  of  Xj)  (density  of  flow  leaving  Xj)P^  Y 
where  X,Y  are  any  combination  of  C,  V  except  the  combination  V,V. 


The  flow  which  leaves  a  spherical  zone  V-  in  the  time  interval 
(n+l)T<T<(n+*2)T  is 

2  ^^2 

2,r, 


(5.4.3) 


Similarly,  the  flow  which  leaves  the  spherical  zone  C.  during  the  time 
interval  (nfl)T<T<(n+2)T  is 


2N 

2«/co8ejp(Cj)^^Z^  - 


2iir^^cosej^p(Cj)^P^  ^  +2*r^^co8ej  (5.4.4) 

*J  X  M 


I 


40 


AEC/C-TDR.63-88 


where 

cosGj 

K  -  1 
=  2 


-  cosCSj-t^) 

r  KKNj^ 

I  Nj^<I<2N 


(5.4.5) 

(5.4.6) 


2  2  2 

Noting  that  r  /r  *  sin  7»  Py  ^ 

V  C  VjOj 

and  denoting  cos6jp(Xj)^  by  a  (^)jj»  Eqs. 
as  follows: 


2 

=  Pp  /sin  7  (see  Equation  5.2.6) 
(5.4.3)  and  (5,4.4)  are  rewritten 


2  2N 

(l-o<  )/sin  7  L  cose  p(C  )  P 

J=1  ^j''l 


(5.4.7) 


(l-a„ 


2N 

S  cos6 
J=1 


(5.4.8) 


where  K  is  defined  by  (5.4.6).  The  flow  chart  for  the  molecular 
transfer  is  self-explanatory.  The  molecular  incidence  for  each 
spherical  zone  is  found  by  the  relation 


2  r/*  2«(Xj)^j 

_ n^O _ 

2 

2  r^^  Z:!ecose^(l-a^  ) 


00 


n”0 

(l-a 


)cosej 


(5.4.9) 


The  infinite  series,  however,  is  very  rapidly  coverging  and  in  most 
cases  it  is  sufficient  to  take  n^lO  or  less. 


5.5  MODIFIED  APPROACH  TO  THE  COMPUTER  SOLUTION 

This  approach  is  formulated  in  matrix  notation  and  has  the  following 
advantages  over  the  approach  described  previously. 

a)  There  are  no  limitations  to  the  number  of  surface  discontinuities 
of  the  parallel  type. 

b)  It  is  sufficient  to  know  only  the  initial  flows.  Thus,  calculations 
of  flows  which  leave  the  surfaces  in  subsequent  time  intervals  are 
unnecessary. 

c)  The  final  solution  to  the  flux  distribution  is  not  the  result  of  a 
chopped  off  summation  of  an  infinite  series. 

d)  It  is  free  of  errors  introduced  in  the  calculations  of  the  flows 
which  leave  the  surfaces  during  nT<t<(nihl)T ,  n>2. 
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where  0^  *  0^+  “ 

being  the  lower  latitude  of  the  I'th  spherical  zone. 
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5.5«2  MoUculor  Transfer 

The  following  difference  equations  are  easily  derived 

2N 


F  (nfl)  *  (l.a„  )  Z  P 
J 


1=1  CiV/Cj(n) 

2N 


F  (rri-1)  =  (1-a  )  E 

S  S  1=1 


P  F  (n)+P  F  (n) 


In  matrix  notation  the  above  two  equations  take  the  form 
F^(rrfl)  =  V'c/c(^> 


Fc(rH-l)  = 


P-VcFv(-'>+P'cC^c("> 


where  the  primes  denote  transpose  matrices* 

The  boundary  conditions  then  are 

Fv(0)  =  2«r/^  AyEF^ 

Fj,(0)  =  2«rg^^  A^EFg 

From  (5. 5.2, 3)  and  (5. 5.2.4)  in  view  of  (5. 5. 2. 5)  and  (5.! 
obtain 

F„(i)  . 


■  *c 


2,r„2ieF'vcVV2”c^"‘"cC*C®fc 


Substituting  (5. 5.2.3)  into  (5. 5.2.4)  we  obtain 

r‘"vcVcv^c<'‘>-"’"cc'c<'«-«' 


Fc(iH"2) 


Summing  up  (5. 5.2.9)  from  n  =  0  to  N  «  oo  and  solving  for  H- 


I-A^ 


P' 


cvV 


CV 


-1 


sin^7 


+  P* 


CC 


F^dH 


V'cvVc/c«» 


,  2 
sin  7 


where  I  Is  the  unit  matrix. 


(5.5.2. 1) 

(5. 5. 2. 2) 

(5. 5. 2. 3) 

(5. 5. 2. 4) 

(5. 5.2. 5) 

(5. 5. 2. 6) 
>.2.6)  we 

(5. 5. 2. 7) 

(5. 5. 2. 8) 

(5. 5.2.9) 

,  we  obtain 
(5.5.2.10) 
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Similarly,  summing  up  (5.3.2.3)  from  n  »  0  to  n  »  oo  we  obtain 

(5.5.2.11) 


“v- Vov  |“. 


5.5.3  Molecular  Flux  Distribution 

2 

Since  the  area  of  the  I'th  spherical  zone  on  surface  X  is  cosG^^ 

where  6^  =  =  6^.+  ^ 


the  rate  of  incidence  on  this  spherical  zone  is 


V 


"Sc,  = 


(5. 5. 3.1) 


Hence,  in  matrix  form 

Then  from  (5.5.2.10)  and  (5. 5. 3. 2)  we  obtain 


CV 


sin^7 


+  P' 


CC 


-1 


where 


p* 


p*  + 
CC 


cA"' 


EA„ 


CV 


,  2 
sin  7 


Similarly,  from  (5.5.2.11)  and  (5.5.3. 2)  we  obtain 
,-l 


sin^7 


V’cv  [ 


EA^Rc+AcEFc 


(5. 5.3.2) 


T  (5. 5. 3. 3) 


(5. 5.3.4) 


,-l 


,  2 
sin  7 


"cv“c 
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6.0  COMPUTER  RESULTS 


For  studying  general  trends  in  the  molecular  flux  distribution  on  the 
chamber- vehicle  surfaces,  a  variety  of  cases  were  chosen  and  the  results  were 
plotted  on  semi log  paper  (see  figures  16-26).  Each  of  these  figures  contains 
two  graphs  for  the  molecular  flux  distributions,  one  for  the  chamber  wall 
(dotted  curve)  and  the  other  for  the  vehicle  (circle  curve).  The  cases  under 
study  were  chosen  in  terms  of  the  following  parameters. 

a)  Size  of  sun  (large  and  small) 

b)  Size  of  vehicle  (large  and  small) 

c)  Pumping  coefficients 

d)  Source  of  outgassing 

The  size  of  the  sun  is  determined  by  the  discontinuous  line  N.  while  the 
size  of  the  vehicle  is  determined  by  the  magnitude  of  N2.  For  the  size  of  the 
sun  and  the  vehicle  there  were  chosen  two  values  for  and 

«  30  Large  size  sun 

33  Small  size  sun 

=  6  Large  size  vehicle  (7=30°) 

=  3  Small  size  vehicle  (7=15®) 

The  different  combinations  of  a),  b),  c),  and  d)  making  up  the  different  test 
cases  are  shown  in  the  following  table. 

TABLE  OF  TEST  CASES 

2N  -  36 


Case 

h 

^2 

a 

V 

h 

q 

m 

1 

30 

6 

.9 

( 

) 

( 

) 

0 

1 

0 

10 

2 

33 

6 

.9 

0 

1 

0 

10 

3 

33 

3 

.9 

0 

1 

0 

10 

4 

30 

6 

.95 

0 

1 

0 

10 

5 

30 

6 

.8 

0 

1 

0 

15 

6 

30 

6 

.9 

1 

0 

0 

10 

7 

30 

6 

.9 

0 

0 

1 

10 

8 

30 

6 

.8 

0 

0 

1 

15 

9 

33 

3 

.9 

0 

0 

1 

10 

10 

30 

6 

.9 

.1 

.1 

.3 

.5 

1 

10 

11 

30 

6 

.9 

.1 

.1 

.01 

.1 

1 

10 
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Looking  at  these  graphs  we  can  point  out  a  number  of  interesting  phenomena. 
Generally  speaking,  in  the  cases  of  a  non-outgassin*;  vehicle  the  flux  distri¬ 
butions  on  the  chamber  wall  and  on  the  vehicle  cross  each  other  at  a  point. 

It  is  not  as  much  the  size  of  the  sun  as  it  is  the  size  of  the  vehicle  which 
affects  this  crossing  point.  A  comparison  between  cases  1)  and  2)  {Figures 
16  and  17)  shows  that  the  crossing  point  has  been  affected  by  .less  than  5 
while  between  cases  2)  and  3)  (Figures  18  and  19)  the  crossing  point  has  been 
displaced  by  20^.  Also  a  slight  increase  in  the  pumping  coefficient  on  the 
chamber  wall  will  hardly  affect  the  crossing  point  as  it  can  be  seen  from 
cases  1)  and  4)  (Figures  16  and  19).  An  interesting  phenomenon  is  presented 
in  case  6  (Figure  21).  The  two  flux  distributions  almost  coincide.  If  the 
area  of  the  sun  shrinks  to  zero  we  have  the  uniform  case  of  out gas sing  chamber 
wall.  Then  the  molecular  flixx  on  the  unpumping  areas  is  both 

the  chamber  wall  and  the  vehicle  as  it  can  be  seen  from  equation  3.2.5. 

Setting  F^=  1  and  *  .9  we  find 


I 

9 


Thus  in  case  6  the  flux  distributions  on  both  surfaces  of  the  chamber  wall  and 
of  the  vehicle  are  equal  and  10  times  as  much  the  flux  distributions  in  the 
uniform  case.  From  this  observation  it  appears  that  the  presence  of  a  non¬ 
pumping  non-outgassing  sun  will  have  a  very  slight  effect  on  the  uniformity 
of  molecular  flux  distribution.  Of  course,  this  cannot  be  the  case  for  all 
sizes  of  vehicle.  From  the  ratio  of  the  area  of  the  sun  to  the  area  of  the 
vehicle  in  case  6  we  can  state  safely  that  both  fliuces  are  equal  to  a  constant 
depending  on  the  size  of  the  sun  as  long  as 


sun  area  ^  2 
vehicle  area  — * 

As  we  have  implied  above,  in  the  cases  of  predominantly  out gas sing  vehicle 
the  two  fluxes  do  not  cross  each  other.  Further,  the  flux  distribution  on  the 
chamber  wall  remains  practically  uniformly  distributed,  apart  from  a  slight 
negligible  variation,  while  the  flux  on  the  vehicle,  increases  from  a  minimum 
at  6  ■  -it/2  to  a  point  between  0^  and  30^  above  the  equation  at  a  constant 
rate  and  from  that  point  on  at  a  higher  rate  leveling  off  around  9  «  n/2. 

Case  7  is  the  superposition  of  cases  1,  6  and  10  each  multiplied  by  an 
appropriate  factor  if  we  neglect  the  pumping  coefficient  and  a^.  Thus  if 
we  multiply  case  1  by  .3  and  case  6  by  .3  and  substract  the^sum  of^these  two 
results  from  case  10  we  will  obtain  approximately  case  7.  To  see  this  let  us 
select  three  points  6  «  -](/2,  0,  jc/2  1+d^  denotes  the  flux  of  the  I'th  case 
at  any  6  we  must  have  d-^j^Q-(.5dj^+.3dg)  ^  either  on  the  chamber  wall  or  on 
the  vehicle.  Consider  che^zlin  on  the  chamber  wall  at  9  •  x/2« 
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-  4.9248x10"^,  -  1.0238,  dj^^  -  5.4276xl0‘^  d^  -  2.4138x10 

dio  -  (.5d^+.3dg)  -  .23316 


At  0  -  0 

d^  -  9.1821x10"^,  dg  -  1.0154,  d^^  -  6.2723x10' 

d^  -  2.9613x10’^ 

dio  -  (.5d^+.3dg)  =  2.7670x10"^ 


At  0  ■  «/2 

d^  -  1.9230x10’^,  dg  -  1.0107,  d^^  -  6.527xlO'^ 

d^  =  2.8306x10'^ 

dio  -  (.5dj^+.3dg)  =  .25334 

The  agreement  is  better  at  d  »  n/2  and  becomes  less  accurate  with  Increasing 
d.  This  is  of  course  due  to  the  fact  that  we  have  discarded  the  pumping 
coefficient  a^2  case  10.  The  flux  distribution  on  the  vehicle  can 

be  considered  in  a  similar  way.  Thus,  if  N-,  a^^,  a^-  and  a  are  kept  fixed, 
the  results  of  any  combination  of  source  outgassing  can  be  obtained  from  a 
linear  combination  of  the  following  three  cases 


i)  =  1,  92  =  9  “  0 

ii)  ■  0,  =  1,  9  =  0 

ill)  q^^  »  q2  ■  0,  q  -  1 

The  superposition  principle  becomes  also  useful  in  the  cases  where  many  types 
of  gases  are  considered  with  independent  pumping  and  outgassing  characteristics 
for  each  type  of  gas.  The  solution  then  in  that  case  is  obtained  by  superim¬ 
posing  the  individual  solutions  for  each  type  of  gas. 
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FIGURE  2 
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FIGURE  6 
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FIGURE  10 
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FIGURE  n 
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FIGURE  12 
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FIGURE  13 
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FIGURE  14 
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CASE  NO.  2  TOTAL  HITS  PER  UNIT  AREA 
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CASE  NO.  4  TOTAL  HITS  PER  UNIT  AREA 
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FIGURE  19 
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CASE  NO.  7  TOTAL  HITS  PER  UNIT  AREA 


FIGURE  22 
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CASE  NO.  8  TOTAL  HITS  PER  UNIT  AREA 


ZONE 


FIGURE  23 
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CASE  NO.  9  TOTAL  HITS  PER  UNIT  AREA 
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FIGURE  25 
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CASE  NO.  II  TOTAL  HITS  PER  UNIT  AREA 


FIGURE  26 
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PROBABILITY  FLOW  CHART 

FROM  CHAMBER  WALL  TO  VEHICLE  AND  FROM  VEHICLE  TO  CHAMBER  WALL 


A  ENTRY 


FIGURE  27 
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PR08ABILITY  FLOW  CHART 

FROM  CHAMBER  WALL  TO  CHAMBER  WALL  WHEN  TT/Z  -Zy  S  0 


eo  TO  0 


FIGURE  28 
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MOIAItLITY  FLOW  CHART 

FROM  CNAMRCR  WALL  TO  CHAHKR  «ALl  WHEN  TT/t  ~ly  >  9 


FIGURE  29 
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MOLECULAR  TRANSFER  FLOW  CHART 


FIGURE  30 
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raOBABILITY  FLOW  CHART  (REVISED) 
FROM  CHAMBER  WALL  TO  VEHICLE 


FIGURE  31 
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PROBABILITY  FLOW  CHART  (REVISCO) 

FROM  CHAMBER  WALL  TO  CHAMBER  WALL  WHEN  y  >  TT/^ 
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FIGURE  32 
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